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We discuss the connection among three distinct classes of models often used to explain the late cosmic 
acceleration: decaying cosmological term, bulk viscous pressure, and nonlinear fluids. We focus on models that 
are equivalent at zeroth order, in the sense they lead to the same solutions for the evolution of the scale factor. 
More specifically, we show explicit examples where this equivalence is manifest, which include some well know 
models belonging to each class, such as a power law A-term, a model with constant viscosity, and the Modified 
Chaplygin Gas. We also obtain new analytic solutions for some of these models, including a new Ansatz for the 
cosmic term. 



I. INTRODUCTION 

The combination of a large set of astrophysical observa- 
tions shows that the Universe is currently undergoing a phase 
of accelerated expansion. A natural explanation for this phe- 
nomenon, in the framework of Einstein's theory of gravitation, 
would be the presence of a cosmological constant. However, 
if this constant is interpreted as the vacuum energy, a mis- 
match of about 50 to 100 orders of magnitude occurs between 
the different contributions to the vacuum energy and the ob- 
served value of A (this is known as the cosmological constant 
problem). Besides, if the cosmological term is constant, it im- 
plies that we live in a very special period of cosmic history 
where the contribution of A to the total energy density is of 
the same order of magnitude of that of the matter fields (this 
is the so-called coincidence problem). 

An interesting alternative, put forward by M. P. Bronstein 
already in 1933 [1], is the possibility of a decaying cosmo- 
logical term. Since then, several phenomenological models 
were proposed for this component. Some authors argue that 
a dynamical A is even a requirement of Quantum Field The- 
ory and provide first principles estimations of its dependence 
on cosmic expansion (see, e.g., refs. |2, \M)- A suitable A 
decaying model could potentially connect the primordial in- 
flation with the present acceleration and also alleviate the co- 
incidence problem. 

Another possibility as a driving force for the accelerated ex- 
pansion is the presence of an isotropic viscous pressure (see, 
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e.g., refs. , which can arise from dissipation processes 

or as a consequence of particle creation. 

A third road to cosmic acceleration is provided by so-called 
exotic fluids, which have negative relativistic pressure, such 
as the the widely used linear equation of state (EOS) p = 
wp, with w < —1/3. However, a fluid with nonlinear EOS 
may have more suitable properties, such as being stable with 
respect to pressure perturbations. An example is given by the 
Chaplygin gas 0, Ufa. A particularly interesting property of 
this type of fluid is the possibility of unifying dark energy 
(DE) and dark matter (DM) (for a review of this scenario, 
see, e.g., ref. |9J]). 

Many other frameworks which lead to an acceleration 
of the cosmic expansion have been proposed, most no- 
tably those based on scalar fields — genetically denoted as 
quintessence. For recent reviews on DE, with special empha- 
sis on quintessence, see, e.g., refs. ITol [Til [l2tl . It is worth 
noting that all these approaches assume that the constant con- 
tribution of vacuum energy plus cosmological constant are 
identically zero due to some compensation mechanism and 
therefore do not address the cosmological constant problem. 

In this work we discuss the connections among the first 
three possible explanations for the accelerated expansion 
mentioned above, namely, decaying A, bulk viscous pressure, 
and exotic fluids with nonlinear EOS. We show that these 
three classes of models modify in the same way the evolu- 
tion of the expansion rate. Therefore for a given model in one 
class, it is possible to find equivalent models belonging to the 
other classes that produce the same expansion history of the 
Universe. 

We discuss some explicit examples of this connection, dis- 
playing the equivalent models for two popular sets of A de- 
caying cosmologies: one where the cosmic term has a power 
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law dependence on the scale factor and another where it has a 
quadratic form in terms of the expansion rate. 

In the process of investigating the equivalent models, we 
derive analytical solutions for the scale factor in some partic- 
ular cases. We also introduce and find analytical solutions for 
a new model of decaying A term. These solutions are useful 
to test numerical codes, they provide some insight for more 
generic numerical solutions, and are handy for pedagogical 
purposes. 

A corollary of the discussion presented in this work is that 
one cannot observationally distinguish these three possible 
explanations for the accelerated expansion using only back- 
ground results (i.e., the zeroth order homogeneous evolution 
of the Universe), such as the supernovae distance-redshift re- 
lation. 

The paper is organized as follows. In section|II]we present 
a brief review of the three frameworks to explain the acceler- 
ated expansion on which we focus in this work, making appar- 
ent the connections among them. Some specific examples are 
discussed in section [Till where equivalent models are explic- 
itly shown and analytic solutions are found. We summarize 
our results in section [IV] Finally, we present our concluding 
remarks, and discuss avenues of future research in section Pvl 

Throughout the text we use natural units, where c = G = 
h = k,B = !■ 



II. THREE POSSIBILITIES FOR COSMIC 
ACCELERATION 

In this section we briefly review the three roads to cosmic 
acceleration and present some models that will be discussed 
throughout this work. We introduce the basic equations that 
govern the evolution of the scale factor in these frameworks, 
exhibiting the formal analogy between them. 



A. Decaying cosmological term 

A dynamic cosmological term (i.e., a component with 
energy-momentum tensor T„ v = A (x) g^) emerges natu- 
rally in the framework of quantum field theory (see, e.g., 
ref. H and refs. therein). Cosmological models with time- 
varying A were introduced in the 1980's (see, e.g., refs. 
lEH El El d ES EH) either as alternatives to the infla- 
tionary paradigm and to a primordial singularity, or to recon- 
cile inflation with observational data. Such A decaying cos- 
mologies attracted more interest in the 1990's, also in con- 
nection with the age problem, and more models where pro- 
posed and tested against the available data (see, e.g., refs. 
lEH S3 EH IH HI Hi El)- Finally, this type of model 
became popular in the turn of the century, impelled by su- 
pernovae and cosmic microwave background (CMB) evidence 
for cosmic acceleration and a smoothly distributed DE. 

In most A decaying models, the density associated to the 
cosmological term is either determined by the dynamics of a 
scalar field or is given as an explicit expression in terms of the 
scale factor a and/or the expansion factor H. Here we shall 



consider only the later case. For example, Ozer and Taha fl4ll 
introduced a model in which A oc a~ 2 . This same dependence 
was proposed by Chen and Wu lfl9ll . based on dimensional 
arguments in the context of quantum cosmology (see also B201 
12611 ). Gasperini, using the thermal interpretation of A 12711 . 
introduced the more general power law form [ 17] 

A oca"" 1 . (1) 

Several observational consequences of this model were inves- 
tigated by Waga and collaborators 1I2TI1 and it was found to be 
consistent with the existing data for m > 1.6. More recently, 
models in which the energy density associated to the cosmic 
term has the form 

epmo _3+ e 
Pa = Pad + 7, a 

3 — e 

where also proposed JHEMl and compared with observational 
data j29l[30n . In particular, a comparison of Large Scale Struc- 
ture (LSS) data with CMB observations leads to the constraint 

e < 0.7 x icr 3 1I29I 

Models in which A has a scaling in terms of the expansion 
factor were also proposed. For example, investigations using 
the renormalization group a ppro ach lead to a running of the 
cosmic term of the form 1 l3i l3ll[32ll 

A = A + vH 2 . (2) 

This model is consistent with a number of observational data, 
but again a comparison with LSS and CMB data implies v < 
2.3 x 1CT 3 [29] (see also JMl)- On the other hand, an estimate 
of the cosmological term from the trace anomaly of quantum 
chromodynamics (QCD) yields [33] 

A = aH , (3) 

where the proportionality constant would be related to the 
QCD cutoff energy scale (er = Aq CD ). 

Several other models, with different phenomenological 
scalings in H, a, or a combination of the two, have also been 
introduced (see, for instance, refs. 2 JH El El El El)- 
However, in what refers to specific models, in the remaining 
of the paper we shall restrict to the simple scalings given by 
eqs. O, ©, and ©. 

Now, let us state the basic equations that govern the dynam- 
ics of models with a cosmic fluid and a cosmological term (in 
a Friedmann-Lemaitre-Robertson- Walker Universe). This 
will allow to make the equivalence with the other two frame- 
works explicit. Since a component with p\ = —pA can only 
be dynamic if it interacts with another matter-energy compo- 
nent, we shall impose only that the total energy-momentum is 
conserved, leading to the energy conservation equation 

p + 3H (p + p) = -p A , (4) 



A quadratic scaling was also introduced previously, motivated by dimen- 
sional arguments 12211 . 
2 See also ref. (341 . for a collection of several phenomenological A-decay 
laws and some historical remarks. 
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where H — a/a and the dot denotes derivative with respect to 
the cosmic time t. Differentiating the Friedmann equation 



o 87T . . k 

H 2 = — (p + pa) - -j 



and using equation leads to 



k 

H - — = -4tt (p + p) 



(5) 



(6) 



From now on, we shall consider only the case in which the 
cosmic fluid has a linear barotropic equation of state: 



P= (7- !)/>• 



(7) 



Using eq. (0 and inserting the expression above into eq. © 
one gets 



H 



47T7/9A ■ 



(8) 



Once an expression for pa in terms of H and/or a is given, 
this equation can be solved to obtain the time behavior of the 
scale factor. 



B. Viscous pressure 

The investigation of possible roles of dissipative processes 
in the Universe has accompanied several developments of cos- 
mology in the past decades. The consideration of models 
with dissipation started to raise considerable interest in the 
1970's, both as a candidate to explain the high entropy per 
baryon ratio inferred from the CMB (see, e.g., ref. IB6IP . as 
well as a mechanism for isotropization and homogenization 
of the Universe (see, e.g., |H3]). Later, viscosity was invoked 
as a way to avoid the primordial singularity (see, e.g., ref. 
1 38]) and to drive an inflationary expansion (see, e.g., refs. 
|3g @(1 H~0, HI]). Its role in the transition from a de Sitter 
epoch to the Friedmann epoch (deflationary phase) was also 
investigated (see e.g. B43l 14411 ). More recently, models with 
bulk viscous pressure were investigated as possible sources 
for the current phase of accelerated expansion of the Universe 
(see e.g. refs. J4j,|5[]). 

In the context of a homogeneous and isotropic Universe, 
dissipation can only be present through a bulk viscous pres- 
sure (also known as second viscosity), whose effect in the 
energy-momentum tensor (and therefore, on the equations of 
motion) is to add a new term to the isotropic dynamic pressure 



P = p-n, 



(9) 



where P is the equilibrium (thermostatic) pressure and II is a 
correction term present in dissipative situations. This term can 
appear either due to a real bulk viscosity — to be derived from 
first principles from kinetic theory (see, e.g., ref. lf36L l45ll ) 
— or as an effective pressure associated to the phenomenon 
of particle creation, which arises naturally in the context of 
quantum processes (see e.g. j46Tl ). 



For first order deviations from equilibrium it ma y be shown 
that the generic form of the extra term is given by I36ll47ll 



n = (9. 



(10) 



where 9 is the divergence of the four-velocitiy field u M , which 
in a homogeneous Universe, for comoving observers, is sim- 
ply given by 6 = 2>H. The viscosity coefficient ( can be a 
function of the dynamic variables (p, p, etc.), but not of its 
(space-time) derivatives i36l l47ll . The second law of thermo- 
dynamics requires that £ > IT36T1 . On phenomenological 
grounds, it is usually assumed that ( is a function of the en- 
ergy density p only. For example, a generic power law form 
for the viscosity coefficient 



C = ap" 



(ID 



was introduced and investigated in [48] and several analytic 
solutions were discussed in j40ll and in subsequent papers. 

As mentioned above, a bulk viscous pressure can also ap- 
pear as a consequence of quantum effects, such as particle pro- 
duction and trace anomaly (for a calculation from first princi- 
ples in the case of a scalar field, see ref. [46]). An Ansatz to 
represent the effects of particle creation, which includes non- 
linear terms in the expansion factor, was given by Novello and 
Araujo 114911 . In their work, the viscous term takes the form, 



n 



JV 

E 

k=l 



(12) 



where the coefficients can be functions of p. 

Prigogini and collaborators 115011 applied the thermodynam- 
ics of open systems to cosmology and derived the bulk vis- 
cosity from particle creation. Assuming that the created parti- 
cles are in thermal equilibrium with the existing ones and that 
the creation process occurs at constant specific entropy, they 
found 3 



n 



P 



n6 



(13) 



where n is the particle number density and ^ is the parti- 
cle source (if ^ > 0) or sink (if \& < 0), which is given 
by W — iV M :M , where A M = nu^ is the particle flux vec- 
tor. Calvao et al. I5II1 give a thermodynamical description 
of particle creation processes in the Universe, still restricting 
to adiabatic transformations, but relaxing the assumption of 
constant specific entropy. They propose the Ansatz 



n = 



(3^ 



(14) 



which generalizes expression ( fT3l . 

As in the previous sections, we shall consider only a fluid 
with the equation of state 4 P = (7 — 1) p. Inserting expres- 



3 This form of Prigogini et al.'s result 1 50] is taken from ref. IBUl . 

4 In models with 7 > 2/3, the bulk viscous pressure can drive the acceler- 
ated expansion if n > (7 — 2/3) p. However, in this case the correction 
term in (9) dominates over the thermodynamic pressure and is not a small 
correction to the energy-momentum tensor. Therefore this term is more 
naturally interpreted as arising from particle creation processes. 
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sion (0 in equation © and using the Friedmann equation (0 
with no cosmological term we obtain 



H 



2 



30 



- 1 = 4?rn . 



(15) 



A comparison with equation I© makes evident that the vis- 
cous pressure can be equivalent to a cosmological term. If 
both II and p\ are expressed as functions of the scale factor 
or the Hubble parameter alone, the explicit equivalence of the 
two models is straightforward. For example, a (linear) bulk 
viscosity of the form ( fTOb with constant viscosity coefficient 
is equivalent to a cosmological term proportional to the expan- 
sion factor: p\ = (3(/j)H. Therefore, Schiitzhold's model 
(eq. 0) for A is equivalent to a linear viscosity with constant 
coefficient. However, in general, if n and p\ have different 
explicit expressions in terms of p, H, or a, the equivalent 
models can only be found after the solution of the equation 
of motion is obtained. 



C. Nonlinear equations of state 

Models with nonlinear equations of state have raised a sub- 
stantial interest recently, mainly due to the possibility of DE 
and DM unification lf52l l53ll . For instance, if p is a nonlin- 
ear function of the energy density p, it could be vanishing in 
high density regions, behaving as DM, and be close to — p 
in low density ones, acting as DE. This behavior would lead 
naturally to a transition from matter domination to a de Sitter 
state in the cosmic evolution and could have a minor impact 
in processes taking place in the early Universe. Moreover, 
such equation of state can lead to a positive (adiabatic) sound 
velocity (c^ = dp/ dp), even for negative pressures. 

An example of an equation of state exhibi ting these proper- 
ties is given by an inverse power law J^HH |53fl 



M 4(qH 



V 



p" 



(16) 



where M > has dimension of mass and a is a real dimen- 
sionless constant. The case a = 1 can be derived from the 
Born-Infeld action for a scalar field [54] and arises from the 
embedding of theories in higher dimensional space-times (in 
particular from (3 + 1)— branes immersed in a (4 + 1)— bulk), 
being connected with string theory H% 1551 l56Tl . Because of its 
similarity with an EOS proposed long ago in the context of 
aerodynamics (8[], this fluid became know as Chaplygin gas, 
and the generic power law ( fT6l ) as Generalized Chaplygin Gas 
(GCG 0, l53ll ). This type of matter component was dubbed 
as quartessence l57ll . since in this framework, only a fourth 
component, besides radiation, baryons, and neutrinos, would 
be needed to describe the cosmic evolution. 

The GCG model is consistent with a number of observa- 
tional data involving the background evolution of the Universe 
for -1/2 < a < 1/2 (see, e.g., refs. HI, HI). On the other 
hand, equation ( fT6l ) fails to reproduce LSS and CMB data if it 
is applied to the fluctuations Q59D . However, if the sound ve- 
locity vanishes (for example, due to the introduction of a cer- 



tain type of intrinsic entropy perturbations), then the model is 
again in agreement with the data [60, 6JJ] . 

Several other nonlinear EOS have been discussed in the lit- 
erature (see, for instance, refs. (§, HI |H HI H HI |M]). 
Effective equations of state appear generically in scalar field 
models with noncanonical purely kinetic Lagrangians 16711 . 
For example, a generalization of the Born-Infeld action II68I1 
leads to equation ( TTBT l. while other generalizations yield new 
classes of quartessence models (see, e.g., ref. 16911 ). 

A natural generalization of equation ( fT6l ), in the context of 
the linear E OS d7 ]) considered in sections dll At and ( 1IIBI ) is 
given by j7Sl2ll|72l|7l|7l 



P= (7- 1)P 



M 4 (" +1 > 



(17) 



which became know as the Modified Chaplygin Gas (MCG) 
EOS. 

This type of EOS appears in several settings with differ- 
ent physical motivations. For example, the consideration of 
Schwarzschild-Anti de Sitter black holes in 5D B75I1 leads to 
the EOS 



2 12 



pP ' 



(18) 



where I is associated to the curvature scale of the asymptotic 
anti-de Sitter space. 

Another particular case of equation ( fT71 ) that is useful in the 
cosmological setting, is given by 



(19) 



which provides an excellent approximation for a relativistic 
gas of massive particles lf76Tl . 

For the purpose of illustrating the equivalence of the three 
frameworks, it is convenient to consider more generic EOS 
with a linear and a nonlinear part 



P= (7- 



(20) 



Inserting this expression into equation (0 and using the Fried- 
mann equation (f5]l with no cosmological term we obtain 



37 



4tt/ (p) 



(21) 



Clearly, the similarity of the equation above with equations 
(O and (15[ is manifest, showing that a variable cosmic 
term, a bulk viscous pressure, and a nonlinear term in 
the EOS can play an equivalent role in the determination 
of the background cosmological dynamics. Although the 
aforementioned equivalence between the three frameworks 
holds for generic equations of state, we have restricted to 
fluids with a linear term to simplify the discussion of some 
specific examples that shall be developed in the next section. 

It is worth noticing how the basic equations are modified 
in these three classes of models. The introduction of a new 
matter-energy component — the dynamical A term — alters 
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both the energy conservation and the Friedmann equations 
(eqs. |4] and [5] respectively) with respect to a single fluid. On 
the other hand, the inclusion of a viscous pressure modifies 
only the energy conservation, through the inclusion of EL Fi- 
nally, a nonlinear term in the EOS does not change any equa- 
tion formally (although it appears in both equations). 



HI. EXAMPLES AND ANALYTIC SOLUTIONS 

In this section we shall illustrate the connection among the 
three roads to cosmic acceleration discussed above through 
some examples, showing the explicit equivalence among spe- 
cific models in each framework. In particular, we shall con- 
sider two common types of dynamical A models: a power law 
cosmic term of the form (Q3 and A-decaying models with a 
quadratic dependence on the expansion rate, including mod- 
els (|2]i and (01 as particular cases. We explicitly find their 
analogues in terms of a bulk viscous pressure and a nonlinear 
EOS in some particular cases. In this process we derive novel 
analytical solutions to these models. 



A. Cosmological term depending on the scale factor 

For models in which the A term is given as a function of the 
scale factor, the equivalence with a model with nonlinear EOS 
of the form ( f2Qb can be easily tested if the energy conservation 
equation allows an analytic solution, providing an explicit ex- 
pression for p (a). In this case, p\ (a) is directly linked to the 
non-linear term / (p (a)). 

As an example, let us consider the following EOS: 



P = (7- 



Bp 



1 + Vi + Cp 

B 



= (7-l)p+^(l-VT+CpJ, (22) 

where the parameter B is dimensionless, and C > has di- 
mensions of inverse energy density. The solution of the energy 
conservation equation @, for 7 7^ 0, 5 is 



Vl + Cpo + l-B/7/ — I — 



VVl+Cpo-1 J 



B-2-f 



-37 



(23) 



Choosing B = 7 this expression can be easily inverted, 
leading to 



(24) 



where a = C' 1 (V 1 + Cp - l) a 37 and c 2 = 2^[/C 
Therefore, the nonlinear term is given by 



/ (p) = % (VT+Cp~ - 1 



-37/2 



' 2 



(25) 

We thus see that this nonlinear term plays the same role as 
a dynamical cosmic term of the form (Q~|l, with m = 3j/2. 
Hence, in this case, the power law A decaying model with a 7- 
fluid (eq. [7]) is equivalent to a single fluid with EOS ( l22l i. with 
B = 7. From equation d25l ), it is clear that these models are 
also analogous to a fluid with linear EOS and viscous pressure 
of the form IT = Hqo,-^/ 2 (see eq. [BJl. 

If on the one hand the evolution of the scale factor will be 
the same in the three models, on the other hand analytical so- 
lutions may be easier to find with one specific choice of the 
equivalent descriptions. In particular, while it is not possible 
to find a simple analytic solution for the energy conservation 
equation (01) with a power law decaying A term and a linear 
fluid, such solution is easily found for the nonlinear EOS ( f22b . 
The advantage of one description over the others in deriving 
the evolution of the system will become even more apparent 
in the discussion bellow. 

Analytic solutions for ait) for the equivalent models dis- 
cussed above can be easily found for 7 = 2rt/3 (where n is 
an integer number). For example, let us consider the "radia- 
tion era" (7 = 4/3), such that in the decaying A model we 
have pa a~ 2 (or, in the case of viscous models, II oc a~ 2 ). 
In this case, equations ([8]), ( fT3T >, and d2"TT l allow the simple an- 
alytic solution 6 



167TCT 



fc r + At 



1/2 



(26) 



where a is the constant of proportionality present in each 
model. The constant of integration A is left free, since equa- 
tions ([8j, ( TOl ). and d2Tb are second order differential equations 
(and we have already used a(0) = 0). 

The solution above is obtained in a simpler way for the non- 
linear fluid model, since in this case we have an explicit ex- 
pression for the energy density as a function of the scale factor 
(eq. 1241 . which can be inserted in the Friedmann equation ©, 
without a cosmological term, yielding 



ff 2 = ^/ 37 +C2a -37M_ A. 



87TC5 



For 7 = 4/3 this equation becomes 



H 2 



87TC1 _ 4 

— — —a 



k 



(27) 



(28) 



which is analogous to the one obtained for a Universe with 
only radiation and curvature, with the quantity k — 871x2/3 in 



' The particular case 7 = leads to the relation 

(VI + Cp - 1) exp VI + Cp 
(VI + Cp - 1) exp VI + Cp ~ 



"0 



3S/2 



6 Here and throughout we shall impose the condition a(0) = 0. On the other 
hand, to simplify the expressions, the "normalization" of an = a(to) will 
be left arbitrary. 
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the place of the curvature parameter k. Now we have a first 
order differential equation which can be easily solved to get 



8 " C2 k]t 2 + ^ l87TCl 



1/2 



(29) 



Now the constant A has a clear interpretation in terms of the 
EOS parameters and the "normalization" p(ao) — po. 

A simple analytic solution can also be found from eq. d27| i 
in the case 7 = 2/3: 



27TC9 o 

-t 2 + t 



8nci 



(30) 



Of course, this solution is also valid for a cosmological term 
decaying as a -1 in a Universe dominated by a fluid with EOS 
p = —p/3 (e.g., a gas of nonrelativistic cosmic strings) for 
any curvature k. 

These simple examples illustrate how the use of equivalent 
models is handy for obtaining analytic solutions. In particu- 
lar, for decaying A models where p\ is a function of the scale 
factor a, it may be more advantageous to use an equivalent 
model in terms of a nonlinear fluid, since it may allow to ob- 
tain analytic solutions for p(a). 



B. Cosmological term depending on the expansion rate 

We shall now consider the case where the cosmic term is 
expressed as a function of H. More specifically, we introduce 
the Ansatz 



A = A + aH + vH 2 



(31) 



which can be regarded as a series expansion of a generic 
A(H). 

Clearly, from equations (0 and ([TBI l. it is easy to see that 
the dynamics of a Universe with a cosmological term of the 
form above and filled by a perfect fluid is equivalent to the 
one of a Universe with no cosmic term, but with and imper- 
fect fluid with viscous pressure of the form (fT2l with constant 
coefficients (and including a constant term). 

On the other hand, to make the comparison with nonlinear 
fluid models — for generic curvature — the combination of 
the Friedmann and energy conservation equations has to be 
used. Some examples that admit analytic solutions for a(t), 
are discussed below. 



1. Trace anomaly cosmological term 



Let us first consider Schiitzhold's model [33] of decaying A 
(eq. [3), which has the form (l3TT l with Ao = and v = 0. As 
mentioned in section ITTBl in the case of a linear fluid (eq. |7J, 
this model is equivalent to a viscous fluid with constant vis- 
cosity coefficient. 

In a flat Universe, this model leads to the analytic solution 
EHl77l|7i 



[exp (<77*/2) - 1] 



(32) 



As pointed out by S. Carneiro 17911 . this same solution is 
found in the case of the MCG (eq. [T7J with a = -1/2. In 
fact the energy conservation equation allows a simple analytic 
solution for this type of fluid 7 : 



P = Po 



a -3 7 (l+«) + 



(33) 



, l+Q 



where po is a constant of integration and A = (M 4 /po) 
Thus, in the flat case and for a = —1/2, the nonlinear term is 
given by 



/ ( p ) = A/y 72 = HM 




(34) 



showing the equivalence of the two models. 

For some values of 7, this nonlinear fluid model with a = 
— 1/2 allows analytic solutions for any value of the curvature 
k. For example, if 7 = 4/3 (which can be regarded as a 
"generalized radiation fluid") one has, 



k ( smhJ^M 2 t 



3M 2 



1 y J^-m 2 



(35) 

while for 7 = 2/3 ("generalized cosmic string gas") one has 

2 

2tt P 1 /2 M 2 1 - 




2P /2 



These new solutions generalize equation (l32l . for these two 
choices of 7, for arbitrary curvature. Although these solu- 
tions are not directly related to a constant bulk viscosity or 
cosmic term proportional to H, they exemplify how the use of 
this nonlinear fluid model, by analytically solving the energy 
conservation equation, simplifies the obtention of analytic so- 
lutions for a(t). 

Analytic solutions for a = 1 and 7 = 4/3, in the flat 
case (i.e., for EOS [T9i>, were obtained in reference |[76ll . It 
is worth mentioning that this model is equivalent (again in the 
flat case) to a linear fluid with particle creation, leading to a 
bulk viscous term of the form (TT~4b with constant source and 



(3 oc p 



-1/2 



On the other hand, if (3 = const., this particle 



creation model is equivalent to a MCG with a = 1/2. 



7 For 7 = the solution is 



P = PO 



1 + 3 1 + a 

\ Po J 



(1+c) 



111 - 
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2. Renormalization group cosmological term 

Another interesting case is found for a = in equation 
d3~Tl i, which gives the form (0 representing the effect of quan- 
tum corrections to the vacuum energy in the renormaliza- 
tion group approach fH HH HH], In the fiat case, the effect 
of this cosmic term is clearly analogous to having a term 
f(p) — (Ao/87r) + (■yi//3)p in equation (OTl i. This analogy 
is still valid for generic curvature, as can be seen by inserting 
expression © into equation ([8]): 



H 



ta 
2 



(37) 



This equation has the same solution as for a model with 
constant A and linear EOS of the form (0, but with 7' = 

7(1 - v/3), k' = fe(37/2 - l)/(7(3 - v)/2 - 1), and 
A' = A /(l-u/3). 

The case 7 = 1/ (3 - 
lytic solution given by 



(i.e., 7' = 1/3), allows an ana- 



a(t) 



sinh 



Ao 
3-v 



T-(2f-3) 



(38) 



An explicit solution for 7 = 4/(3 — v) can also be obtained 
as a particular case of equation d44l . presented in the next 
section. 



5. Quadratic expansion rate cosmological term 

As a final example, let us consider the full quadratic form 
(eq. [3Tb . Inserting this expression in equation ([H) and intro- 
ducing the new variable u, such that a = u n , with 



(39) 



7(3-1^) 



one has 



.. 7^ • 7 A o 

u — u — u 

2 2n 



±(^1-1)^ ,40, 



For k = (or 7 = 2/3), the solution is 

.— P 70-*/4 



it = uh 



+c 2 exp 




, (41) 



where c\ and C2 are arbitrary constants. 

Imposing a(0) = as in the previous cases, we have 



a(t) 



sinh 



7 2 a 2 
16 



7A0 
2n 



(42) 



Notice that for v = Ap = we recover the solution (1321 . as 
it should be. Also, for a = and 7 = 1/(3 — v), we recover 
expression (l38l for fc = 0. 

Solutions for k ^ can be found in some particular cases. 
For example, for n = 1/2, such that 7 = 4/(3 — v), the 
solution is 



,37 ,\ 2 

U = 1iff + K — -1 — j— 

2 / 7A0 



Again, imposing a(0) = 0, we have 
a(i)= je^*/ 4 sinh 
2k /37 



(43) 



, 7^ 
* [ ~~r~z + 7A0 



1G 



7 A C 



■-1 l-e 7CT * /4 cosh 



>(44) 



A nonlinear fluid model that leads to the same solutions 
as the quadratic A(H) model, in some cases, is given by the 
following EOS 



p=( 7 -l)p-M 2 p 1 / 2 -A., 



(45) 



with A being a constant. This model allows an analytic solu- 
tion for the energy conservation equation, which, for 7 / 0, 
is given by 8 



Si 



+ 1 



-3 7 



1/2 

Pa 



m 2 n 

2T V i 



(46) 



l M 2 . fl 



/ J () 



1/2_ M 2 1 / A/4 T A 
2 T i^" 1 " 7 



Clearly, for fc = the model above (eq. |45I > is equivalent to 
the quadratic A model (eq. |3T1 >. and has a solution of the form 
(l42l . However, this solution would hardly be found from in- 
verting relation d46*i i and inserting the result in the Friedmann 
equation. This is an example where the use of a nonlinear 
EOS makes it harder to obtain analytic solutions, contrary to 
the examples of sections lIII Al and flHB II 



IV. SUMMARY 

In this work it was shown that three distinct classes of cos- 
mological models (decaying A, bulk viscous pressure and flu- 
ids with exotic nonlinear equations of state) are equivalent, in 
the sense that in each of the classes there are models that re- 
produce exactly the same expansion history of the Universe. 



' For 7 = the relation is 

A 



aV2 



1/2 

■Pa 



M 2 



In 



M 2 p 1 / 2 + A 
M 2 pl /2 + A 



3M 2 



In — . 

ao 
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This equivalence is apparent from the formal analogy among 
equations ©, ( fT51 ), and d2~TT). 9 

As explicit examples where the equivalence is manifest, we 
have considered two popular classes of "decaying vacuum" 
cosmologies: a power law A term (sec. IIII Al l and a cosmic 
term depending on the expansion rate (sec. IIIIBl i. 

We have analyzed the case of a generic quadratic cosmolog- 
ical term of the form A = Ao+aH+vH 2 , which is equivalent 
to a viscous term of the form ( TTSl i with N = 2 and constant 
coefficients, plus a constant term. In the flat case (and for a 
linear EOS p = (7 — l)p) this model is also equivalent to the 
MCG Gas with a = — 1/2 plus a constant term (eq. |45l >. 

We obtained a complete analytical solution for these mod- 
els in the flat case (eq. 14Tb . For arbitrary curvature, we have 
shown a solution for 7 — 4/(3 — v) (eq. l44l . As far as we 
know, it is the first time that such type of model is discussed 
in the literature. 

For Ao = 0, and in the flat case, this model is equivalent to 
a linear viscosity of the form (fTOb with coefficient C = Co + 
dp 1 / 2 and, for a linear EOS, to the MCG (eq. |T7j with a = 
-1/2. 

For Ao = = v we recover the "trace anomaly cosmo- 
logical term", which is equivalent to a linear viscous term 
with constant viscosity coefficient. Again, in the flat case and 
for a linear EOS, this model is equivalent to the MCG with 

a =-1/2. 

We have found two analytic solutions for the MCG with 
a = —1/2 and arbitrary curvature: one for the "generalized 
radiation fluid" (eq. l35l l and another for the "generalized cos- 
mic string gas" (eq.l36ll. This generalizes the solution for the 
flat case discussed in [78]. 

For (j = we recover the "renormalization group cosmic 
term". We have show that this A-decaying model is equivalent 
to a model with constant A and a linear fluid, with redefined 
values for the cosmological constant, the curvature, and the 
equation of state parameter (see sec. IIII B 2b . We have ob- 
tained an analytic solution for the case 7 = 1/(3 — v) and 
arbitrary curvature. 

Regarding the power law A term, we found that mod- 
els with a linear equation of state p = (7 — l)p and a 
power law cosmic term of the form A oc a~ 37 / 2 (or equiv- 
alently a viscous pressure of the form II = IIoa~ 37 / 2 ) are 
equivalent to a model with EOS given by p = (7— l)p + 
7C _1 (l — VI + Cp) and no A. We have displayed the ana- 
lytical solutions for a (t) for 7 = 4/3 and 7 = 2/3, for generic 
curvature. 

The use of equivalent models can be handy in finding an- 
alytic solutions. In fact, in some cases one particular model 
excels the other ones in the task of deriving explicit solutions. 



9 It is worth mentioning that the equivalence among these frameworks has 
been discussed in previous works for some particular cases. For example, 
the equivalence between variable-A cosmological models and models with 
bulk viscosity was pointed out in refs. I8(JI (for one specific model in 
the flat case) and | 81] (for two models). Also, the similarity of the trace 
anomaly cosmic term with the MCG with a = — 1/2, in the flat case, was 
shown in ref. (7^1 . 



Here, we have illustrated this property with particular exam- 
ples where the solutions to a particular model are easily found 
using the equivalent model belonging to another class. 



V. DISCUSSION AND CONCLUDING REMARKS 

As mentioned above, a practical application of the mathe- 
matical equivalence pointed out in this work is in the search 
for analytical solutions, which can be useful for pedagogical 
purposes and to test numerical codes. But, besides this oper- 
ational use, does the equivalence among the three classes of 
models have a deeper physical connection? Can we observa- 
tionally distinguish these models? 

Certainly they cannot be distinguished with observables 
that are mostly dependent on the background cosmology, such 
as the redshift-distance diagram of type la supernovae. On the 
other hand, although their dynamic role is equivalent, these 
models may have different thermodynamic al behavior. For 
example, it is know that bulk viscosity and matter creation 
processes are thermodynamically distinct (see, e.g. B82IP 10 . 
Thus, one may try to use thermodynamic arguments to distin- 
guish among the three frameworks. 

The existence of equivalent models in the three classes dis- 
cussed in this work is not very surprising given the limited de- 
grees of freedom of a homogeneous Universe. In fact a simi- 
lar equivalence was pointed out in refs. 118-411 connecting scalar 
field models with nonlinear fluids in the homogeneous case. 

A way to distinguish these models is to study their behavior 
when more degrees of freedom are present. Conversely, if 
the equivalence among the three frameworks remains in more 
generic configurations, this could point to a deeper connection 
among them. Thus, the natural next step in this investigation 
is to study the behavior of linear perturbations in the three 
classes of models. 

Naturally, if A decaying models are considered homoge- 
neous by construction, the perturbations of the inhomoge- 
neous component will be clearly different from viscous mod- 
els or models with nonlinear equations of state, since the vis- 
cous or nonlinear terms are intrinsically inhomogeneous. On 
the other hand, the vacuum contribution can be regarded be- 
ing as space-time dependent. For example, in models where 
the cosmic term depends on the expansion rate, H can be re- 
placed by a third of the local value of the velocity divergence 
(see, e.g., H). 

Models with bulk viscous pressure are readily applicable to 
inhomogeneous situations by using the local values of the den- 
sity and velocity divergence in the expression of the viscous 
term. Of course, models which are identical in the flat case for 
homogeneous configurations will have a different behavior in 
the perturbations. 



The thermodynamics of a vacuum decaying model is discussed by Alcaniz 
& Lima 1 30], whereas the thermodynamic behavior of the GCG was inves- 
tigated by Santos, Bedran, & Soares t§3l . 
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It is also straightforward to consider adiabatic perturbations 
in models with nonlinear equations of state. However, for 
quartessence models such as the GCG, this type of perturba- 
tions generate large scale structures that are inconsistent with 
observational data (see, e.g, B59L 16311 ). On the other hand, a 
specific type of entropy perturbation brings these models back 
in agreement with the data 1 60, [6jl [fiH] . Several other ways to 
generalize to the inhomogeneous case equations of state that 
correspond to the GCG for the background have also been 
considered (see, e.g., refs. SUllMl). 

The study of perturbations in the three frameworks is rather 
involving and presents several subtleties. Therefore, this anal- 
ysis is left for subsequent works (see, e.g., B89I0 . 

It is also worth investigating the microphysical motivation 
of the three frameworks. For example, it is known that models 
with nonlinear equations of state may arise from scalar fields 
with noncanonical kinetic terms lf67L I68L [69h . 

Another avenue of research is to investigate wether the 
equivalence holds with respect to other frameworks for cosmic 
acceleration, such as self-interacting gases 0, Ht] and other 
dark energy models. 

Finally, several cosmological observables (such as large- 



scale structure, cosmic microwave background fluctuations, 
and distance-redshift relations) have to be obtained in order to 
check the viability of these models with respect to the diverse 
set of observational data available. 

Once these steps are undertaken, the study of equivalent 
models may prove useful for finding physical insights and mo- 
tivations to phenomenological and first principle models that 
provide a good description of the available astro-cosmological 
data. 
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